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Abstract 

The gravitino may well play an important role in cosmology, not only because 
its interactions are Planck-suppressed and therefore long-lived, but also because 
it is copiously produced via various processes such as particle scatterings in ther- 
mal plasma, and (pseudo) modulus and inflaton decays. We study a possibility 
that the early Universe was gravitino-rich from various aspects. In particular, a 
viable cosmology is possible, if high-scale supersymmetry is realized in nature as 
suggested by the recent discovery of the standard-model like Higgs boson of mass 
about 125- 126 GeV. We find that the Universe can be even gravitino-dominated, in 
which case there will be an entropy dilution by the gravitino decay. If the gravitino 
abundance is proportional to the reheating temperature, both the maximal baryon 
asymmetry in leptogenesis and the dark matter from the gravitino decay become in- 
dependent of the reheating temperature. The dark matter candidate is the Wino-like 
neutralino, whose mass is suppressed compared to the anomaly-mediation relation. 
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1 Introduction 



The recent discovery of the standard- model (SM) like Higgs boson of mass about 125 - 
126 GeV [U[2] may imply high-scale supersymmetry (SUSY) [31 H]. SUSY relates fermions 
to bosons and vice versa, and adds the superpartner of each particle in nature; hence, in 
a supersymmetric theory there will be a fermion that partners with the graviton, called 
the gravitino, which becomes massive via the super-Higgs mechanism after SUSY gets 
spontaneously broken. The interactions of the gravitino are suppressed by either the 
Planck scale or the SUSY breaking scale, and therefore it is naturally long-lived. 

The gravitinos are produced in various processes. In the early Universe, the gravitinos 
are produced by particle scatterings in thermal plasma [5j El [7J. Also, it is known that 
any scalar field with a non-vanishing vacuum expectation value (VEV) generically decays 
into gravitinos, if kinematically allowed. In particular, a large amount of gravitinos can 
be produced by the (pseudo) modulus [H [91 [10] and inflaton [HI [TJl [131 CHI EES] decays. 

If the gravitinos are abundant in the early Universe, they may play an important role 
in cosmology, because of their longevity. Indeed, it is well-known that too many gravitinos 
significantly affect the evolution of the Universe in contradiction with observations; the 
cosmological gravitino problem [HI El HE]- For instance, if the gravitino mass is around 
the weak scale, it typically decays during big bang nucleosynthesis (BBN), altering the 
light element abundances in contradiction with observations. In the case of high-scale 
SUSY, however, the gravitino can decay before BBN, and the success of the standard 
BBN remains intact. Thus, the high-scale SUSY may lead to a viable cosmology even if 
many gravitinos are produced in the early Universe. The main purpose of this paper is to 
study such gravitino-rich Universe from both cosmological and phenomenological aspects. 

What is peculiar to the gravitino is that it is coupled to any sectors in nature. Thus, 
once produced, the gravitinos will distribute their energy to all the lighter particles in- 
cluding the SM particles. For instance, the gravitino decay produces the lightest ordinary 
supersymmetric particle. If it is the lightest supersymmetric particle (LSP) and if the 
R-parity is conserved, the LSPs thus produced may account for the observed dark matter 
(DM) abundance. Alternatively, if there are light degrees of freedom in a hidden sec- 
tor, the gravitino necessarily decays into them. Some of the hidden-sector particles thus 
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produced may contribute to DM, if they are long-lived and non-relativistic. If they are 
relativistic at the BBN and recombination epoch, they will increase the expansion rate, 
which may account for the dark radiation hinted by recent observations [22| [23]. (See e.g. 
Refs. pi E51 ES E3 EEJ ESI M EI] for the models.) Also, if the inflation took place in 
a hidden sector, the only way to reheat the visible sector may be through the gravitino 
production; the gravitino may play a role of messenger because of its universal couplings. 
Thus, the gravitino-rich Universe seems to have various interesting implications. 

In this paper we shall study if the Universe becomes gravitino-rich by considering 
several production processes of the gravitinos; thermal production as well as non-thermal 
production by the (pseudo) modulus and inflaton decays. We find that the gravitino- 
rich Universe is realized for a wide range of parameters, and the Universe even becomes 
gravitino-dominated in some cases. In the latter case, there is an entropy dilution by the 
gravitino decay. Interestingly, the gravitino abundance and therefore its dilution factor are 
proportional to the reheating temperature in the case of thermal production and the non- 
thermal production by the pseudo-modulus decay. Then, the maximal baryon asymmetry 
in thermal leptogenesis becomes independent of the reheating temperature. Assuming 
the Wino-like LSP of mass several hundred GeV, the LSP abundance can account for 
the observed DM. For this, the anomaly mediation relation for the gaugino masses must 
be modified, requiring a mild cancellation of order 10% in the Wino mass. Note that 
both the baryon asymmetry and the neutralino LSP abundance are independent of the 
reheating temperature in this case. 

The rest of this paper is organized as follows. In Sec. El we briefly discuss the gravitino 
lifetime and abundance. We study various production processes of the gravitinos in Sec. El 
In Sec. H] we discuss cosmological aspects of the gravitino-rich Universe, and study the 
maximal baryon asymmetry in the thermal leptogenesis as well as the DM abundance. 
In Sec. we discuss how the gaugino mass relation can be deviated from the anomaly 
mediation relation. The last section is devoted for discussion and conclusions. 
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2 Gravitino abundance 



Let us summarize here the abundance and lifetime of the gravitino. For the moment 
we assume the minimal particle content, namely, the minimal supersymmetric standard 
model (MSSM) particles and the gravitino. Later we shall consider a case that there are 
light degrees of freedom in the hidden sector. 

The decay rate of the gravitino into the MSSM particles is approximately given by 

r 3/2 «, s/2 - MSSM) K ig^, (1) 

where m 3 / 2 is the gravitino mass, Mp ~ 2.4 x 10 18 GeV is the reduced Planck mass, and 
we have assumed that the gravitino is heavier than all the MSSM particles, for simplicity. 
The decay temperature of the gravitino, T 3 / 2 , is defined by 

T 3/2 , (^)^, (2) 



« „. UGeV (£^)- (_£*_)•, (3) 

where g*(T) counts the relativistic degrees of freedom at temperature T. If the gravitino 
is heavier than several tens TeV, its decay temperature is above a few MeV, avoiding the 
tight BBN constraint [19]. We adopt m 3 / 2 = 10 3 TeV as a reference value in the following. 

Let us quantify the gravitino abundance in terms of the ratio of the gravitino number 
density to the entropy density; 

^ = =sa, (4) 

where Sj represents the entropy density originated from the inflaton decay. Note that 
does not include the entropy produced by the gravitino decay. 

If the gravitinos are abundant, they may come to dominate the energy density of the 
Universe, since the gravitino is massive and long-lived. The gravitino-dominated Universe 
is realized if 

> eo 

4m 3 / 2 
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This is equivalent to the condition that the entropy produced by the gravitino, s 3 / 2 , be 
larger than the pre-existing entropy, Using Eq. (j3J), this condition can be rewritten as 



^ 2 - L5xl ° l^^J lionw) ■ (6) 

If this condition is met, the entropy dilution factor A is given by 

A _ s 3/2 + Sl ^ n 3/2 f n ^Y\ (7) 

Si Si \ S3/2 / 

~ 90( 9 * {T3/2) Y ( m3/2 Y*(M\ ( 8 ) 

where s 3 / 2 represents the entropy density produced by the gravitino decay, and we have 
used (713/2/53/2) = 3T 3 / 2 /4m3/2- Noting that Y3/2 ~ 10 -2 when gravitinos are in thermal 
equilibrium, one can see that the dilution factor varies from 1 to ~ 10 5 for the gravitino 
mass around 10 3 TeV. 

It is possible that a modulus field dominates the Universe after the reheating, and 
decays into gravitinos. If the branching fraction of the modulus decay into gravitinos is 
large enough, both the modulus and gravitino decays produce entropy. In this case the 
total entropy dilution factor is the same as Eq. (JS}, and it can be also expressed in terms 
of the modulus abundance as 



Si PX \S3/2 

% 2 N 3/2 ^(^\, (10) 



'm x Ty 2 V s i 

where S3/2 denotes the branching fraction into gravitinos, iV 3 /2 is the average number 
of gravitinos for the decay modes into gravitinos, mx and px are the mass and energy 
density of the modulus X. When the modulus mainly decays into a pair of gravitinos, we 
expect B 3 / 2 — 1 and N 3 / 2 = 2. 

If ([6]) is not satisfied, the gravitinos do not dominate the energy density of the Universe, 
but they can still have a significant impact on cosmology. For instance, if its abundance 
is greater than ~ 10~ 12 , the LSPs produced by the gravitino decay tend to overdose the 
Universe, as long as the annihilation of the LSPs is not effective and the LSP mass is larger 
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than or comparable to (9(100) GeV. We therefore adopt the criterion for the gravitino-rich 
Universe as, 



Gravitino-rich Universe : Y3/2 > 10 12 



(11) 



Of course, there are various possibilities, for instance, the R-parity may be violated, and 
the QCD axion may be the dominant component of DM, or there may be other lighter 
SUSY particles in the hidden sector, and so on. A different criterion should be adopted 
depending on the scenarios of interest. To be concrete, however, we adopt the above 
criterion in this paper. 

3 Production processes 

Now we consider various processes for the gravitino production. 
3.1 Thermal production 

Gravitinos are produced by particle scatterings in thermal plasma, and its abundance is 
given by [3 El [7] 

where we have omitted the logarithmic dependence on Tr as well as terms that depend 
on the gaugino masses. The definition of Tr is given by 

t rS p^Tpyf^, (is) 

where denotes the total decay rate of the inflaton, and we assume that the infla- 
ton mainly decays into the MSSM sector. Thus, the gravitino abundance increases in 
proportion to the reheating temperature. 

The gravitino-rich Universe is realized if the reheating temperature satisfies 

T R > 10 10 GeV, (14) 

while the gravitino dominates the Universe if the reheating temperature is 

rH > 10 , Ge v(^)-*(^. m 
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Thus, an extremely high reheating temperature is required for the gravitino to dominate 
the Universe. 

Note however that the expression for the gravitino abundance (fl2l) does not take ac- 
count of the longitudinal mode. The contribution from the longitudinal mode, namely the 
goldstino, becomes relevant if the gauginos are heavier than the gravitino. Furthermore, if 
the inflaton decays into the SUSY breaking sectoral the would-be goldstino may be ther- 
malized, leading to Y$J^ ~ 10~ 2 , where we have assumed that the relativistic degrees of 
freedom is of order 10 2 . Thus, it should be kept in mind that the above estimate on the 
gravitino abundance ({12]) depends on the assumption about the reheating process and the 
mass spectrum. 

3.2 Non-thermal production 

Now let us consider non-thermal production of gravitinos. In the following we consider a 
pseudo modulus in the SUSY breaking sector, a general modulus field, and an inflaton as 
possible source of gravitinos. 

3.2.1 Pseudo modulus field in the SUSY breaking sector 

We first consider a pseudo-modulus decay into gravitinos. A pseudo modulus in the 
SUSY breaking sector is a plausible candidate for the gravitino production because of the 
following reasons; (i) the existence of the pseudo modulus is generic in the SUSY breaking 
models; (ii) the coherent oscillations are produced because the potential minimum for the 
pseudo modulus during inflation is generically deviated from the low-energy minimum; 
(iii) it predominantly decays into a pair of gravitinos. 

We study the low-energy effective theory of O'Raifeartaigh type SUSY breaking model. 
After integrating out the massive fields, the Kahler and super-potentials are 



where z is a pseudo modulus field, A is a cut-off scale, n represents the SUSY breaking 



W 



K 




(16) 
(17) 



scale, and the constant Wo — m^^Mp is fixed so that the cosmological constant almost 



1 The inflaton itself may be a part of the SUSY breaking sector. 
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vanishes at present. We can assign an R-charge 2 on z, which is explicitly broken by Wo 
down to Z 2 r- All the parameters are set to be real by an appropriate redefinition of z 
and the U(l)# transformation. 

The z is stabilized at the potential minimum, 

(*> = (18) 

where m z is the mass of z given by 

m z = (19) 

The F-term of z is given by F z ~ — /i 2 ~ v / 3^3/2^p, and SUSY is indeed broken. The 
precise value of A depends on details of the SUSY breaking. We simply assume here that 
A is much smaller than the Planck scale so that z is heavier than the gravitino. 

Suppose that the U(l)# symmetry remains a good symmetry during inflation and 
that the inflation scale, H^, is larger than m z . Then, if z acquires a positive Hubble- 
induced mass, z is stabilized near the origin during inflation, and starts to oscillate with 
an amplitude of (z) when H ~ m z after inflation^ Thus the pseudo-modulus abundance 
is given by 

e*~3T R (^X, (20) 



where we have assumed that the pseudo modulus starts to oscillate before the reheating. 
On the other hand, if the Hubble parameter during inflation is smaller than m z , the 
pseudo-modulus abundance gets suppressed. If the inflaton at the potential minimum is 
heavier than m z , the pseudo modulus cannot follow the change of the potential minimum, 
and the coherent oscillations are induced [321 [33] . The abundance is given by 

* * StJ^V^V- (21) 
Si \m z J \m z J 

If the inflaton mass at the potential minimum is lighter than m z , no coherent oscillations 

are induced. In the following analysis we use the pseudo- modulus abundance ( 1201) . As 

one can see from Fig. HJ this is the case for the inflation scale larger than 10 9 ~ 10 10 GeV, 

which covers many inflation models. 



2 Even if the Hubble parameter during inflation exceeds the dynamical scale, it is possible that z is 
stabilized near the origin, when the Hubble parameter becomes so small after inflation that the low-energy 
description ([TC]) and (TIT)) are valid. 
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Figure 1: In the shaded (red and blue) regions, the Universe becomes gravitino-rich (see 
Eq. (Hip ). The gravitino dominates the Universe in the lower shaded (red) region. The 
pseudo modulus also dominates the Universe in the region below the dashed line. We have 
fixed Tr = 10 9 GeV, but the lines are only weakly dependent on Tr. See the text for details. 



The quartic coupling in the Kahler potential induces the z decay into the goldstino 
pair with the decay rate 



V(z ->• zz) ~ ^-^t- (22) 

As long as z does not have any sizable couplings with the SM sectorjf] it predominantly 
decays into a pair of gravitinos. Thus, the gravitino abundance from the pseudo-modulus 
decay is given by 

'(lo^v)(l^v) 4 (To^v)^ <*> 

The Universe becomes gravitino-rich if is greater than 10~ 12 , namely, 

-.SlxlCGe V (J^f (^f. (24) 

3 Planck-suppressed couplings do not change the argument. 
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The gravitino-dominated Universe is realized if 

m z < lxl0 8 GeV (*)* (-^L-) l ° . (25) 

V80/ \WGeX) V10 3 TeV/ V ; 

In Fig. [TJ the Universe becomes gravitino-rich in the shaded (red and blue) regions. In the 
lower shaded (red) region, the gravitino dominates the Universe, while it does not in the 
upper shaded (blue) region. The dotted (purple) lines show the contours of the entropy 
dilution factor, A = 10 3 and 10 6 . In the upper left region, the condition ( |24|) is not 
satisfied and the gravitino abundance is small, while the decay into a pair of gravitinos is 
kinematically forbidden in the lower right region. In the region below the dashed line, the 
pseudo modulus dominates the Universe before the decay. We have fixed T R = 10 9 GeV 
in the figure, but the conditions on m z are only weakly dependent on Tr; both the upper 
and middle solid lines (the dashed line) shift in proportion to oc (T^ 13 ). 

3.2.2 Modulus field 

Next we consider a modulus field as a source of gravitinos. In superstring theories, 
moduli fields necessarily appear at low energies through compactifications. Most of these 
moduli must be stabilized in order to get a sensible low-energy theory, since the moduli 
determine all the physically relevant quantities such as the size of the extra dimensions, 
physical coupling constants, and even the SUSY breaking scale. Many moduli fields are 
stabilized by flux compactifications [311 EH], and the remaining ones can be stabilized a 
laKKLT [20]. 

The detailed properties of the modulus depend on the compact geometry, brane config- 
urations, and stabilization mechanism. Here we consider a modulus field X stabilized a la 
KKLT; it has a mass heavier than the gravitino, and its F-term is suppressed by the mass 
as Fx ~ m 3 / 2 (X) /mx, where m>x denotes the mass of X. During inflation the potential 
of X is considered to be deformedj^ Then the potential minimum may be deviated from 
the low-energy minimum, leading to coherent oscillations after inflation. We assume that 
the energy stored in the coherent oscillations of X dominates the Universe. The entropy 
density Si should be interpreted as that from the modulus decay in this subsection. 

4 There is an upper bound on the inflation scale in order to avoid the destabilization and run-away. 
Even for the inflation satisfying the bound, the cosmological moduli problem still exists [33J ■ 
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We also assume that the modulus X has Planck-suppressed couplings to the visible 
sector so that its decay rate can be expressed as 



* = (26) 



where c is a numerical constant of order unity. The decay temperature of the modulus is 



tt 2 <7*(T 



Tx = Vr X M P , (27) 

- 3G -(^)^(t^v) § - (28) 

It was shown in Refs. [H [9j [10] that such a modulus generically decays into a pair of 
gravitinos and the branching fraction is sizable: B3/2 = (9(0.01 - 0.1)1^ Then the gravitino 
abundance from the modulus decay is estimated as 



Y $ = (29) 



= axio-^d)-^) 1 . (30) 

Thus, the Universe becomes gravitino-rich for a wide range of the modulus mass and the 
branching fraction of the gravitino production. In particular, if the branching fraction 
into gravitinos is of order 0.1, the gravitino even dominates the Universe. (Compare ( 1301) 
with ©.) 

3.2.3 Inflaton 

In a series of works [9j [T0l [TTJ [12j [131 CHI EE [21] , it was revealed that the inflaton generally 
decays into gravitinos. The gravitino production rate depends on whether there is a SUSY 
breaking field that is singlet under any symmetries. Such a singlet is required in the gravity 
mediation to generate gaugino masses of the correct size, while it is not necessary in the 
gauge and anomaly mediation. We assume that there is no such a singlet for the moment. 

The gravitinos are produced by various processes; (a) the gravitino pair production [HI 
[TUl [TTl [12] ; (b) decay into the SUSY breaking sector at tree level [13] ; (c) anomaly- induced 



5 Precisely speaking, this is the case if the SUSY breaking field is heavier than the modulus mass, or 
if there is a SUSY breaking field that is singlet under any symmetries. 
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decay into the SUSY breaking sector at one-loop level [H] . The gravitino production rate 
can be expressed as 



where is the inflaton mass, and ((f)) a VEV of the inflaton. Here it should be noted 
that ((f)) is evaluated at the potential minimum after inflation. The precise value of the 
numerical coefficient x depends on the production processes, possible non-renormalizable 
couplings in the Kahler potential, and the details of the SUSY breaking models [T5]. To be 
concrete, let us assume the minimal Kahler potential and the dynamical SUSY breaking 
(DSB) with a dynamical scale A. As we have seen before, the SUSY breaking field z 
can acquire a mass m z heavier than m 3 / 2 , and we assume m z ~ A ~ y'm^Mp in the 



For a low-scale inflation model with < A, the process (a) becomes effective, and 
x = 1. On the other hand, for the inflaton mass larger than A, the processes (b) and 
(c) become effective instead. The inflaton decays into the hidden quarks in the SUSY 
breaking sector via Yukawa couplings of the hidden particles (process (b)), or into the 
hidden gauge sector via anomalies (process (c)). Since the hidden quarks and gauge 
bosons (and gauginos) are energetic when they are produced, they will form jets and 
produce hidden hadrons through the strong gauge interactions. The gravitinos are likely 
produced by the decays of the hidden hadrons. We denote the averaged number of the 
gravitinos produced per each jet as N^/ 2 - Then x is given by [15] 



where and ah are the Yukawa coupling and a fine structure constant of the hidden 
gauge group, respectively, N y denotes a number of the final states for the process (b), N g 
is a number of the generators of the gauge group, and Tg and T r are the Dynkin indices 
of the adjoint representation and the matter fields in the representation r. Although x 
depends on the structure of the SUSY breaking sector, its typical magnitude is 0(1O~ 2 — 

6 The gravitino production can be suppressed if the SUSY breaking field is lighter than the inflaton 
and the inflaton mass is below the dynamical scale [21] [36] . 
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Figure 2: The gravitino-rich Universe is realized in the region above the blue solid line, 
while the gravitino dominates the Universe in the region above the red solid line. We have 
set x = 1, T R = 10 9 GeV and m 3/2 = 10 3 TeV. 



10) for > a]3 To be concrete we adopt x = 1 as a reference value in the following, but 
it should be kept in mind that there are uncertainties in the gravitino production rate in 
this case. 

The abundance of gravitinos from the inflaton decay is therefore 



Y, 



3/2 



3/2 



3T, 



R 



~ 1 x 10 7 a; 



9*{Tr 
200 



(0) y (_hh_\ 2 

10 15 GeV/ U0 15 GeV/ V 10 9 GeV 



(33) 
(34) 



Thus, the gravitino abundance increases as the inflaton VEV and mass. Also, it should be 
noted that the gravitino abundance is inversely proportional to the reheating temperature. 
The Universe will be gravitino-rich, if 



((f)) > 3 x 10 12 GeV x 



9*(Tr 
200 



m 



4> 



10 15 GeV 



10 9 GeV 



7 Roughly, we expect N 3/2 = 0(1 - 10 2 ), N g = 0(1), a h = 0.1 - 1, and T, 
strongly depends on the SUSY breaking models. 
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(35) 
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and the gravitino dominates the Universe if 

X (lO^Gev) (lO^v) ' (36) 

The regions satisfying (|35|) or ( 136|) are shown in Fig. [2j together with representative 
inflation models, new (371 [38] , hybrid [39] . smooth hybrid [ID] and chaotic |JT] inflation. 
(The effect of the constant term in the superpotential on the inflaton dynamics (121 H3] 
is not taken into account in this figure.) We have set x = 1, Tr — 10 9 GeV and m 3 / 2 = 
10 3 TeV. One can see that the Universe can be gravitino-rich for a large portion of the 
parameter space, and this region will be larger for lower Tr. For instance, in the smooth 
hybrid inflation, the inflaton mass and VEV are of order 10 15 GeV. In the chaotic inflation 
without a Z2 symmetry, the inflaton mass is about 2 x 10 13 GeV, and the VEV is around 
the Planck scale. For such high-scale inflation models, the resultant gravitino abundance 
is so large that the gravitino dominated Universe can be realized. 

So far, we have assumed the absence of the SUSY breaking singlet field z. When 
it is present, it generally enhances the gravitino production rate through the operators 
like K D \<p\ 2 z, or \cf)\ 2 zz. The former induces a kinetic mixing between the inflaton and 
z, while the latter induces the pair goldstino production. Furthermore, such a singlet 
z may be produced as coherent oscillations, and it decays into gravitinos. Thus, the 
gravitino-rich Universe will be more plausible in the presence of such a singlet. 

4 Cosmology 

In this section we discuss cosmology of the gravitino-rich Universe. In particular we study 
the neutralino DM produced by the gravitino decay, and estimate the baryon asymmetry 
in thermal leptogenesis in the presence of entropy dilution. 

4.1 Neutralino dark matter abundance 

We assume that the lightest neutralino \ is the LSP. It decouples from thermal bath at 
the freeze-out temperature, Tf ~ m x /28 for the neutralino mass m x around or above 
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a few hundred GeV. Since T 3 / 2 is lower than Tj for the gravitino at the PeV scale, 
the neutralino abundance produced by gravitino decay is determined according to the 
Boltzmann equation, 



dt 



+ 3Hn x = -{a x v Te[ )n (37) 



where n x is the number density of x- The effective annihilation cross section of the lightest 
neutralino is parameterized by 



cr x v re i) = %r. (38) 



m 2 



One then finds the relic abundance of \ to be [HJ HH 

H 



n x\~ l _ ( n x 



t=t 3/2 + \s(a x v re i) 



(39) 



T=T 3/2 



When the gravitino abundance is large enough, the annihilation among the neutralinos 
will take place until the expansion rate H becomes equal to the annihilation rate. This 
results in 

^ ~ 2 1 x 1(T 12 ( 9 J^A) ^ 1 ( m * ) 2 ( m ^ \ " 3/2 (40) 

s~ V 80 J \10- 2 J v300GevMl0 3 TeVy ' [ ) 

for Y 3/2 > 10- 12 . 

In the gravitino-rich Universe, the lightest neutralino needs to have a large annihilation 
cross section in order not to overdose the Universe. This implies that, even if there are 
neutralinos thermally produced before the gravitino decay, their thermal relic density will 
be small. Let us consider the case where % is dominated by the neutral Winojj] The 
Wino-like neutralino has c x around 10 -2 , where we include the coannihilation effect but 
it is not so large because T 3 / 2 is similar in size to the mass difference of the charged 
and neutral Winos. Fig. |3] shows the parameter region where the Wino-like neutralinos 
from gravitino decay account for the observed DM density. For Y3/2 > 10~ 12 , the neutral 
Wino with mass of several hundred GeV can constitute the DM of the Universe. We 
also plot the contours of the entropy dilution factor A for the gravitino-domination case, 



8 See Ref. [46] for the higgsino DM in the heavy gravitino scenario, and also Ref. [47l [48] for the 
neutralino DM from gravitino decay in the case with Y 3 / 2 less than 10 -12 . 
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Figure 3: Wino DM density from gravitino decay consistent with the observation, for 
m x = 0.3, 0.4, 0.5 TeV, respectively. The left side of the solid line for given m x is excluded 
by the overclosure constraint. The gravitino dominates the Universe in the shaded region. 

Kj/ 2 > 10~ 7 . One can see from the figure that the Wino DM density is independent of the 
gravitino abundance in the gravitino-rich Universe, and that the correct DM abundance 
can be explained for the Wino-like neutralino with mass of several hundred GeV and 
the gravitino at the PeV scale. In the next section we will discuss how to realize such 
hierarchy between the gaugino mass and the gravitino mass. 

Meanwhile, since the charged Wino is nearly degenerated in mass with the neutral 
Wino, Sommerfeld effect should be included when the Wino has mass around Airm\y / g 2 ~ 
2 TeV jl9]. However the resulting annihilation cross section would not be sufficiently large 
to avoid the overclosure constraint for the gravitino mass at the PeV scale. 

4.2 Baryon asymmetry in leptogenesis 

One plausible way to generate the baryon asymmetry is through leptogenesis [50J, and 
the simplest realization is to use lepton number violation of the right-handed neutrinos in 
the seesaw mechanism |51j. Here we consider thermal leptogenesis in the gravitino-rich 
Universe. 
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The observed baryon asymmetry of the Universe can be explained through thermal 
leptogenesis if the reheating temperature after inflation is high enough, Tr > 10 9 GeV 
[52] . On the other hand, the gravitino abundance increases in proportional to Tr when 
produced by thermal scattering or pseudo-modulus decay. Hence one may worry about 
the overproduction of DM by the gravitino decay for the reheating temperature required 
by the successful thermal leptogenesis. In order to avoid the tension, Tr must be chosen 
to be close to about lO 9 " 10 GeV for m 3/2 = 10 TeV- 100 TeV§ 

Such a tension can be relaxed if the neutralino LSP has mass of several hundred 
GeV, while the gravitino mass is at the PeV scale, as the neutralino DM abundance 
becomes independent of the gravitino abundance. In addition, if the gravitino dominates 
the Universe, there is an entropy dilution. When the gravitinos are produced by thermal 
scatterings or the pseudo-modulus decay, the entropy dilution factor is approximately 
proportional to the reheating temperature, implying that the present baryon asymmetry 
reads 

- = 3X10 A ^-jj-^j^—^cc-, (41) 

where k denotes the efficiency factor, m u3 the heaviest neutrino mass, and S e g the effective 
CP phase. Here we consider the lepton asymmetry generated by the decay of the lightest 
right-handed neutrino with mass Mi(< Tr), which is converted to the baryon asymmetry 
via the sphaleron effect. Noting that the maximum value of the baryon asymmetry is 
realized for Mi ~ Tr, and that the dilution factor scales as A cx Tr, one can see that the 
maximum value of baryon asymmetry from leptogenesis will be independent of Tr. 

Thus, if the gravitino- dominant Universe and the mass hierarchy between the LSP and 
the gravitino are realized, both the maximum baryon asymmetry in thermal leptogenesis 
and the neutralino DM abundance are independent of Tr. This feature has an interesting 
implication that thermal leptogenesis can be implemented at high reheating temperatures 
while avoiding overproduction of DM. 

Let us see the case where the gravitino produced by pseudo-modulus decay dominates 
the Universe. In Fig. 2J we show the contours of Wino DM density consistent with the 

9 Interestingly, the neutrino mass anarchy and the observed neutrino mass squared differences suggest 
Tr = 10 9 GeV - 10 10 GeV, when combined with thermal leptogenesis. [53] 
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Figure 4: Gravitino-rich Universe in the pseudo-modulus scenario with m z = 4 x 10 8 GeV. 
The solid lines correspond to the contours of Wino DM density from gravitino decay consis- 
tent with the observation, for m x = 0.3, 0.4, 0.5 TeV, respectively. The gravitino dominates 
the Universe in the shaded region. For successful thermal leptogenesis, the reheating tem- 
perature after inflation should be higher than about 10 9 GeV. The maximum value of the 
baryon asymmetry is constant along the (red) solid line in the shaded region. 



observation in the (777.3/2, Tr) plane, for the pseudo-modulus with mass m z = 4 x 10 s GeV. 
Along the solid (red) lines in the shaded region, the gravitino dominates the Universe, and 
the maximum value of the baryon asymmetry in thermal leptogenesis is constant because 
A oc Tr. For heavier (lighter) m z , the pseudo-modulus abundance decreases, and the 
contour of the entropy dilution factor will shift upwards (downwards). 

5 Mass hierarchy between gaugino and gravitino 

We have seen in the previous section that the Wino-like neutralino LSP of mass several 
hundred GeV produced by the decay of the gravitino at the PeV scale can account for 
the observed DM density. In this section we discuss how to realize such mass hierarchy 
between the Wino-like neutralino and the gravitino. 

So far we have not specified how the SUSY breaking is mediated to the visible sector. 
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In the gravity mediation, both gaugino and scalar masses arise from coupling with a SUSY 
breaking singlet, and all the superparticles have masses comparable to the gravitino mass. 
In this case, the Wino mass can be much lighter than the other SUSY particles only at 
the price of fine-tuning the coupling. Since the naturalness is not a serious issue in the 
high-scale SUSY, we cannot exclude this possibility, and indeed such fine-tuning may be 
required by the anthropic condition on the DM abundance. 

In the absence of such SUSY breaking singlet, on the other hand, the gaugino mass 
can be hierarchically lighter than the gravitino mass. In the heavy gravitino scenario, the 
anomaly mediation gives dominant contributions to the gaugino mass [54"t [55]: 

M A ~ Tg, (42) 

unless one considers a specific SUSY breaking scenario such as no-scale type models. 
Scalar masses are more mo del- dependent; the scalar fields obtain masses ~ ^3/2 through 
contact interactions with a SUSY breaking field in the Kahler potential, which are allowed 
even if the SUSY breaking field is charged under some symmetry. However, if the visible 
sector is dynamically sequestered from the SUSY breaking sector, it is possible that the 
scalar masses are suppressed. Compared to the other soft terms, therefore, the gaugino 
masses are less model-dependent, and have a simple form [56] . 

To arrange a light Wino of mass several hundred GeV for the gravitino mass at PeV 
scale, the anomaly-mediation relation (142 j) is not sufficient, and one needs further sup- 
pression|^j One way is to introduce SM vector-like matter fields \I> + ^ which obtain 
masses from the VEV of a SM singlet so that their loops generate soft masses via gauge 
mediation. Let us consider Ny pairs of ^ + ^ forming 5 + 5 representation of SU(5) with 

W = yySVy + f(S), (43) 

for which the gauge coupling unification of the MSSM is not spoiled. Including the 
contribution from the messenger loops, low energy gaugino masses read 

M a b a F s 

-J a = loV™ 3 / 2 - 16^5"' (44) 



10 In no-scale models, the chiral compensator auxiliary field is much suppressed compared to TO3/2, 
and consequently anomaly-mediated soft masses are much smaller than m 3 / 2 /87r 2 . This provides another 
way to have relatively light gauginos in the heavy gravitino scenario. 
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where (61, 62, ^3) = (33/5, 1, —3) are the MSSM beta function coefficients. 

The anomaly-mediated Wino mass can be cancelled when S is stabilized in such a way 
that (F s I S) /m 3 / 2 is positive and order unity. A simple way to achieve this is to consider 
a composite S having an Affleck-Dine-Seiberg superpotential [57] , 

A 3+n 

f( s ) = -fr> (45) 

where A c is the dynamical scale much lower than the GUT scale, and n is a positive rational 
number. Then the F-term scalar potential from the above superpotential competes with 
the associated A-term to fix S at (A c /m^/2) 1 ^ n+2 'A c with 



F s 2 



rn 



Lx- m 3/ 2l (46) 



S n + 1 V ^ m 3/2 
where m 2 s is the soft mass squared of S. Assuming that S acquires soft mass around 
m 3/2 /8n 2 , one finds that the Bino and Gluino have masses around m 3 /2/87r 2 whereas the 
Wino obtains a relatively light mass for n = 2N^, — 1: 

which is a few hundred GeV for the gravitino at the PeV scale. Since S is a composite, 
the messengers have a small Yukawa coupling yq, ~ A c /Mp, and get a SUSY preserving 
mass M$ ~ (A c /??i3/2) 1 ^ n+2 - ) A 2 /Mp. Thus a dynamical scale higher than about 10 12 GeV 
leads to M# 3> m 3 / 2 as required for the relation (JHEj) to be valid. 

Alternatively, in order to fix 5 at a vacuum giving a positive (F s / S)/m 3 /2, one may 
consider the model with 



K = \Z\ 2 + V( |$| a + ^i^L|$| 



k$ |Z| 

= y£(M 2 -SS'), (48) 

for $ = {E, S, S', ^} and ^> m 3 / 2 , so that S is stabilized near the F-flat direction, 
SS' = M 2 . Here Z is the SUSY breaking field with F-term ~ m 3 / 2 Mp, and k$ = 1 if $ 
resides in a sector sequestered from the SUSY breaking sector. For y^, of order unity, the 
ratio between the soft scalar masses of S and S' is naively estimated by 

mf, (8n 2 ) 2 (l -K S ,)+y 2 

m 2 ~ ( 8 7r 2 ) 2 (l- K5 )+2/|' 1 J 
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Let us assume m| 3> m 2 s , , which would be the case when n s is less than one while k$ = 1 
for the others, or when = 1 for all the involved fields and y <C y^. Then the scalar 
potential develops a minimum at \S'\ 2 ~ (m|/m|,)M 2 3> IS*) 2 with 

?=( 1+ <i)K' (50) 

and therefore the anomaly-mediated Wino mass is cancelled for = 1, making it possible 
to have ~ m 3 / 2 / '(87r 2 ) 2 for m 2 s larger than about 8n 2 m 2 s ,. Meanwhile, the messenger 
mass reads M^, ~ (m|,/m|) 1 / 2 M. 

It is interesting to note that the above model possesses a global U(l), under which 
S and S' have charges of opposite sign while E is neutral. This can be identified as the 
Peccei-Quinn (PQ) symmetry solving the strong CP problem [5E1EH], with the axion scale 
fixed at F a ~ \S'\. The axion also constitutes the DM with Q a ~ 0.4# 2 (F a /10 12 GeV) L18 
where 9i is the initial misalignment. The axion relic energy density will be diluted if the 
gravitino decays after the QCD phase transition, for which F a larger than 10 12 GeV is 
allowed even for large initial misalignment. 

Let us move to the possibility of arranging a light Wino through the higgsino loops. 
Including higgsino threshold effects, the gaugino masses are written 



M a b a k a 

-^ = T ^ms / i + a— 2 m 3/2 , (51) 



for (hi, &2j ^3) = (3/5, 1, 0), and a defined by 



fi sin 2(3 m\ f\\i\ . 

a = — TH~2 ln \ZZ2 » ( 52 ) 



m 3 /2 \fi\*-mx \m; 

where \x is the higgsino mixing parameter, is the mass of the heavy Higgs bosons, 
and tan/3 is the ratio of the Higgs VEVs. For a large \i around 7713/2 and low tan/3, the 
higgsino correction can cancel the anomaly-mediated Wino mass, leading to a light Wino 
at the price of fine-tuning. 

Finally we briefly discuss phenomenological aspects of the model, where the lightest 
neutralino x an d the lightest chargino are both Wino-like. For large fj,, the dominant 
mass splitting between \ and x ± comes from gauge boson loops, and makes the chargino 
heavier than the neutralino by about 160 MeV for around several hundred GeV 
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[6l] . Hence x^ 1 dominantly decays into \ plus soft 7r ± with a decay length of a few cm, 
or longer if boosted, while producing a visible track in the detector. Meanwhile, there 
are cosmological and astrophysical constraints arising since the neutral Wino has a large 
annihilation cross section into a W-boson pair, which require the Wino mass to be larger 
than about 300 GeV [62]. 

6 Discussion and Conclusions 

Being the superpartner of the graviton, the gravitino is coupled to all the sectors with 
Planck-suppressed interactions. Thus, once the gravitino dominates the Universe, it de- 
cays into all the lighter degrees of freedom, including those in the visible sector. Then the 
question is how to generate the baryon asymmetry and DM. In SecJUwe have considered 
thermal leptogenesis and the Wino-like neutralino LSP non-thermally produced by the 
gravitino decay. Here we briefly discuss other possibilities. 

One is the so-called gravitino-induced baryogenesis [63] . In the presence of an R-parity 
and baryon- number violating operator, UiDjDk, and the associated A-term with a CP 
phase, the gravitino decay can generate a right amount of the baryon asymmetry. Since 
the coefficient of the R-parity violating operator must be sizable for the mechanism to 
work, any pre-existing baryon asymmetry would be washed out in this scenario. Also, 
we need a source of CP violation; if there is a SUSY breaking singlet field, it can satisfy 
the requirement. The presence of such SUSY breaking singlet generally enhances the 
gravitino production from the inflaton decay, and furthermore, its coherent oscillations 
can decay into gravitinos. Thus, the presence of the SUSY breaking singlet makes both 
the gravitino dominance and the gravitino-induced baryogenesis plausible. 

Another one is the Affleck-Dine mechanism [01]. The AD mechanism is so efficient 
that a sufficient amount of baryon asymmetry can be generated even in the presence of 
entropy production by the gravitino decay. The AD mechanism in high-scale SUSY was 
explored in detail in Ref. [35] . 

One of the plausible candidate for DM is the QCD axion. If the gravitino decays after 
the QCD phase transition, the QCD axion abundance is diluted, allowing a larger value 
of the axion decay constant [6B] . 
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So far we have assumed that the gravitino is heavier than all the MSSM particles. 
However, it is also possible that the gravitino decay into some of the MSSM particles 
is kinematically forbidden. In this case the gravitino lifetime becomes longer, and the 
gravitino mass is required to be heavier in order to avoid the tight BBN bound. 

If there are additional light degrees of freedom such as the axion and axino, or any 
light degrees of freedom in the hidden sector, the gravitino will decay into those degrees 
of freedom. The branching fraction into these hidden particles can be sizable especially 
if some of the MSSM particles are heavier than the gravitino. Interestingly, there is an 
argument that the string theory contains a plenitude of string axions [66J . If the gravitino 
decays into those axions and axinos with a sizable branching fraction, the produced axions 
(and axinos) can account for the dark radiation hinted by the recent observations [24]. 
Also, if some of the hidden sector particles thus produced are stable in a cosmological 
time scale, they will contribute to the DM; depending on its mass it can be hot or warm 
DM. 

In this paper, we have studied the gravitino cosmology in high-scale SUSY as suggested 
by the recently discovered Higgs boson with mass 125 - 126 GeV. We have discussed var- 
ious gravitino production processes such as thermal production as well as non-thermal 
production by the (pseudo) modulus and inflaton decay, in order to see if the gravitinos 
are abundantly produced in the early Universe. We have shown that the Universe can be 
gravitino-rich for a wide range of parameters, and it is even possible for the gravitino to 
dominate the Universe. In the gravitino-rich Universe, the neutralino LSP must have a 
large annihilation cross section, since otherwise the LSPs produced by the gravitino decay 
would overdose the Universe. We focused on the Wino-like LSP and estimated its abun- 
dance as a function of the gravitino mass and abundance. We also discussed how to realize 
the required mass hierarchy between the lightest neutralino with m x = (9(100) GeV and 
the gravitino at the PeV scale. 

Among various possibilities, we have found that the pseudo-modulus in the SUSY 
breaking sector is a plausible and interesting source of gravitinos. The dilution factor by 
the gravitino decay is proportional to the reheating temperature in this case, which makes 
the maximum value of the baryon asymmetry in thermal leptogenesis independent of the 
reheating temperature. 



23 



If high-scale SUSY is realized in nature, perhaps the naturalness argument is not the 
right guiding principle to understand the observable parameters in our Universe. The 
SUSY breaking scale suggested by the SM-like Higgs boson mass ranges from 10 TeV up 
to PeV or even higher, but we do not know the reason why it takes such a value. It may 
be due to a competition between the bias in the landscape and the anthropic condition 
such as inflation [38] or entropy dilution by the modulus decay [67J. 

The gravitino, being the superpartner of the graviton, is long-lived and has universal 
Planck-suppressed couplings. There are various processes in which gravitinos are pro- 
duced. Therefore the gravitino likely plays an important role in cosmology, especially if 
high-scale SUSY is realized. Our Universe might have been gravitino-rich. 
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